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In This Issue... 

The Primal Heterotope presents a remarkable 
graphic journey, starting from the Mandelbrot set 
and ending at the quasi-Mandelbrot set associated 
with z— 2 + c. On the yellow brick road 
between them, we encounter strange and wonder- 
ful beings, some beautiful, some grotesque. 


Sound and Ceramics reports on Bart Lynch’s 
efforts to generate pleasing pottery from pleasing 
sounds, and llerations on Clay Pots raises the ques- 
tion of extending these ideas to fractal processes. 


The Mandelbrot Set Chaosometer gives a proposed 
measure of the amount of chaos associated with 
the iteration track of a given point in the plane, 
and discusses how well it works as a detector of 
chaos in Mandelbrot calculations. 


1/f Noise discusses a type of fractal noise which has 
proved to be quite useful in modelling certain ran- 
dom processes, and presents an algorithm for gen- 
craling 1/f noise. 


Quicksilver Starship discusses a novel idea for 
drastically decreasing the compute time for gener- 
ating a video or movie of zooming/panning in the 
Mandelbrot set. 

Quaternion Square Root derives a formula for the 
square root of a quaternion. 


The Primal Heterotope 
— Rollo Silver 


A couple of years ago, I started to wonder how to 
deform the quadratic Mandelbrot set (defined by 


ZE 2+ Cc) into the cubic Mandelbrot set 


a 
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(ze 23+ Cc). Could that be done? How? What 
would it look like? 


The topologists call a continuous deformation of 
one figure into another a homotopy; if such a 
deformation exists, the figures are homotopic, and 
topologically equivalent. 


Since the quadratic M-set is connected, while the 
cubic Méset isn't, they aren’t topologically equiva- 
lent, hence they aren’t homotopic. A deformation 
of the one into the other may exist, but it can’t be 
continuous, and cant be a homotopy. I decided to 
call the thing I was looking for a heterotopy. 


My explorations ended up with the one-parameter 
family of Mandelbrot sets M, associated with the 
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| l 2 3 
polynomials P{z) = (1-A zZ +tz +c 
(O<1<1) ,a heterotopy, or non-continuous 
deformation of the quadratic Mandelbrot set M, 


into the cubic set M, . 


I call this particular heterotopy the primal one 
because it is the simplest one I can think of: the 
most primordial one. Simple though it is, it is filled 
with strange and marvelous phenomena. For 
instance, there are certain: singular values of ¢ near 
which the M, change character very rapidly — 
perhaps infinitely rapidly. In hindsight I might 
have expected that. Since the deformation can’t be 
continuous, but ¿ries to be, there must be critical 
points at which the discontinuity expresses itself in 
no uncertain terms. 


In his article 7he Three-dimensional Mandelbrot 
Set (Algorithm issue 3.3, pages 20-21), Earl F. 
Glynn discussed a heterotopy based on the one- 
parameter family of non-polynomial functions 


ze z*+ c, for 1<as3. This approach produces 
an interesting family of images, but is less simple 
than the polynomial approach. For one thing, the 


functions 27+ Cc are multi-valued, which leads to 
peculiar slip-planes in their images. 

In this article, the figures illustrating the journey 
from M, to M, are all greyscale images using one 


of three “color” schemes; they are all 650 x 489 
pixels and have the same escape radius: R= 10. 


The main sequence consists of a series of views all 
having the same magnification. There is contro- 
versy over the term magnification as applied to 
fractal images. Several people have objected to my 
altempt to give a precise meaning to the term: the 
reciprocal of the radius of the largest circle that will 
fit within the field of view. To avoid controversy I 
have coined a new term for this reciprocal radius: 
suidar. Most of the images in the main sequence 
have suidar 0.237 — a value small enough so that 
all significant portions of every image appear on- 
screen. Table 1 (page 9) gives the parameters for 
each of the 31 figures. 

With the time-like parameter ¢ linking together an 
infinity of M-sets, the ensemble of figures achieves 
a sense of motion and change that transcends the 
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static nature of the individual images. 


We start, then, with M, — the familiar Mandelbrot 


set — with t-value 0. Its appearance is familiar — 
except that with the large value of R, the set is 

small when the surrounding regions of dwell 0 and 
1 are shown. 


| Figure 1:ż = 0; suidar = 0.158 


We proceed uneventfully from t= 0 to t= 0.072: 
the dwell boundaries change shape slightly, the 
head of the set grows bigger, and the principal 
buds migrate to the right toward the cusp of the 
cardioid. All is smooth. 


Figure 2:1 = 0.072; suidar = 0.237 
The cubic term does not intrude topologically until 
t becomes comparable to I/R = 0.1. In fact, 
something mildly interesting begins to happen for t 
= 0.073, and is clearly manifest by 0.075: A region 


of dwell 2 breaks through at the critical point -10 
and expands eastward. This region is an artifact of 
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using an escape radius to determine dwells. I do 
not believe that it would appear with pictures 
based on the potential. 


AS 
SES 
= X 


By ¢= 0.08, the intruder has been joined by cohorts 
of dwell 3 and 4, and the normal region of dwell 2 
is elongating westward toward it. 


By t= 0.10 (Figure 3) the cubic term has asserted 
aggressively: the two D = 2 regions have 
merged, a real intruder has appeared, and two 
other lenticulars have appeared north and south. 


=-0.12) 


Figure 3:4 = 0.10 (Suidar = 0.13 


gure 4— a magnified view of the intruder— 


reveals two facts: that the intruder is genuine 


fractal, and that M,, is disconnected: the contour 


younding the region D = 3 is evidently closed, 
enclosing the midget but not the main set. 
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the intruder grown and approaching the main set. 


The Ugliest Fractal: A Digression 

Now let's digress to get a closer look at the 
intruder. Figure 6 shows the intruder with magnifi- 
cation 14 relative to Figure 5. 

I wonder what that funny double structure to the 


left looks like? Figure 7 shows its upper half at rel- 
alive magnification 125. 


He’s my candidate for the ugliest fractal! But like 
the toad, ugly and venomous, he wears yet a pre- 
cious jewel in his head. Figure 8 details the jet 
ejected downward from the cusp between the two 
components of the U.F. at relative magnification 
950. The jet carries an unusual figure: the afore- 
mentioned jewel. 


Figure 8: The Jet 
(t = 0.12; RM = 950) 


Figure 9: The Jewel 
(t = 0.12; RM = 8400) 


Figure 9 shows that the jewel in the center of Fig- 


ure 8 carries a midget in its heart. 


Back To The Main Sequence 


Let’s go back to the main sequence! Continuing 
from Figure 5 Ct = 0.12 ): as t¢ increases the 


intruder approaches the main set— what happens 
as they collide? 


Figure 10 shows the intruder approaching the set 
more closely. A collision is imminent. 


Figure 10: The Intruder Menaces The Set 
(t = 0.1258) 


Figure 11 shows a magnified view (RM = 7.5) of 
the area of imminent collision. 


Figure 11: Imminent Collision 
t = 0.1258; RM = 7.5) 
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The Collision 

This is one of the critical points where things 
change extremely rapidly as fadvances. With / = 
0.12595 (Figure 12) the collision begins. 


Figure 12: The Collision Begins (¢ = 0.12595) 


The collision proceeds as the two Mandelbrot bod- 
ies move together, squeezing the material between 
them into fewer and larger balls (Figure 13). 


Figure 13: Compression and Inflation (t = 0.126 


As the squeeze continues, the balls coalesce into a 
single ball, which becomes the common head for 
the wo Mebodies. Small buds appear at the top 
and bottom of the ball (Figure 14). 


Further compression merges the two bodies and 
the head into a single mass. The buds at top and 
bottom enlarge (Figures 15). Note how rapidly the co Y o ee 
figures change with small changes in f. Figure 16: Separation (t = 0.127) 


Now a surprising thing happens (Figure 16). The 
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lesces the com- 


ed balls forms between the colliders. 
‘compression coa 


y -SIZ 
ponents (Figure 21). 
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achieved. 


Before pushing the rest of the way to t= 1.0, let's 


back off and have a look at the big picture (Figure 
29), 


Figure 28: Almost Complete (¢ = 0.50; RM = 7.5) 


SS 


Figure 26: Ithyphallic Hypertrophy 
(t = 0.32; RM = 7.5) 


Figure 29: Almost Complete: The Big Picture 
(t = 0.50) 


Fi 


gure 27: The Enantiomorphic Urge 
(t = 0.40; RM = 7.5) 


Figure 30:The Cubic Mandelbrot Set, Complete 
(t = 1.00) 
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Figure 30 shows the cubic Mandelbrot set com- 
plete, centered in the field of view. 


And finally, an enlarged view of the Cubic set's 
head: 


Figure 31:The Cubic Set's Head 
(4 = 1.00; RM = 7.5) 


It’s curious that almost all the action takes place in 
the first alt. 087205. the second 


half,0.5 <¢< T is sedate by comparison. 


All the figures were generated using special soft- 
ware written to generate them. 


Table 1: Parameters for Figures 


FIG 
(CLUT) 


suidar center 


0.0 2.0 

0.072 22.0 

0.10 20 

0.10 -7.940 

0.12 230 

0,12 0.131 

0.12 -6.292+0031i 

0.12 -6.2878+ 
0.0064i 

-6.28734+ 
0.00640i 

0.1258 20 

0.1258 195 

().12595 -3.753 

0.126 -3.753 

0.1262 3,753 

0.1266 3.793 


Ore 
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Table 1: Parameters for Figures (Continued) 


center 


-3.753 


0.1285 -3.660+0.548i 


0.1285 -2.0 
0.14 ; -2.0 
0.142 ; -2.0 
0,17 ; -2.0 
O22 -2.0 


0.28 2.0 

0.28 -1.093-1.591i 
0.30 -1.011-1.591i 
0.32 -0.950-1.575i 
0.40 -0.663-1.539i 
0.50 -0.370-1.446i 
0.50 2 

1.00 0 

1.00 178. 115321 


Sound and Ceramics 
Bart Lynch 


Reprinted from LEONARDO Electronic News 


In architecture, natural harmonies occur in Renais- 
sance structures. Harmonic relations of form and 
space were often based on the golden section and 
the ratios therein. These same ratios occur in the 
growth patterns of flowers, fish, and other compo- 
nents of nature. I am currently concerned with 
understanding why these ratios occur and why 
they are pleasing to us. 


I have been translating sounds to three dimen- 
sional pottery using several computer programs in 
order to see if pleasing sounds make pleasing pot- 


tery and vice versa.! Using the sound program 
SOUND EDIT PRO, I can get a visual representa- 
tion of a sound that is time dependent. That visual 
is saved as a picture and imported to SWIVEL 3-D 
where the sound form can be lathed to resemble 
pottery and used as a template to create actual 


1. This work was developed at Deep Creek School in 
conceptual collaboration with Dan Collins. 
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ceramic works. Using these programs, I have also 
been animating the figures so that on the computer 
screen, the pottery forms dialogue with the sounds 
that created them. I see these processes as data 
gathering exercises that help me to understand the 
nature of the harmonic relations so that I will be 
able to use them more effectively in the future. 


Iterations on Clay Pots 
Stewart Dickson 


I have just read Bart Lynch’s piece, Sound and 
Ceramics, in LEONARDO Electronic News 3(7) 
July 15, 1993). This is a very cool piece, on levels 
which go quite a bit deeper than he described in 
his article, I think. 


What he is describing is an iterative procedure in 
which he is translating a physical, 3D object (a clay 
pot, thrown on a potter’s wheel) into an abstract 
representation. He is “sounding” or “ringing” the 
object, recording this event and creating an alter- 
nate, spatial abstract representation of the object: A 
graph, presumably a frequency spectrum) (he did 
not say precisely what graph he produced on 
SOUND EDIT PRO). Wait a minute, on reading 
more closely, he says, “...that is time dependent.” 
This is merely an amplitude envelope, then. A pity. 


At any rate, he then takes the sound envelope 
graph and uses it as a contour template for throw- 
ing a new pot. 


If he were to make frequency spectrum pictures, 
he would be translating the event into the complex 
phase plane and then we could compare the pro- 
cedure to those like computing the Mandelbrot set. 
I think that were Bart Lynch to approach it in this 
way, that he could map out the divergence and 
convergence of his procedure; i.e., he would be 
finding convergent and divergent lathed shapes. 


[ will refer this to some experts in the field (FRAC- 
L) and see what they say. Pm sure the procedure 
could be simulated. I think the acoustical mechan- 
ics of a lathed shape are well known. The question 
is mainly what kind of graphs will produce conver- 
gent sound-producing shapes. I think we can fid- 
dle with the mapping a little until we find 
something that can produce regions of conver- 
gence amid a field of divergent shapes. One can 
imagine that convergent shapes could be analo- 


gous to some kind of purity or “music of the 
spheres”. 


Please let me know your response to these ideas. I 
will compile all responses I receive. 


Stewart Dickson 

The Post Group 

6335 Homewood Avenue 

Los Angeles, CA 90028 

Phone: (213)462-2300 X 832 

FAX: (213)464-1953 

Internet: tpg!dickson@rhythm.com 


The Mandelbrot Set Chaosometer 
—Jay Hill 


Those of you who have followed my e-mail posts 
about the area of the M-set and related computa- 
tions will have noticed the recurrent theme of how 
to avoid long computations by early identification 
of three classes of orbits for the z 22 + c itera- 
tion. These are the escapers, periodics and chaotics. 
The sci.fractals. FAQ (Frequently Asked Questions 
posted on the Internet sci.fractals newsgroup) 
states that the Lyapounov exponent can be used to 
identify chaos. If the Lyapounov exponent A is 
positive, then the iteration is chaotic. The exponent 
is defined as 


Now for the Mandelbrot Set this is 


N 
l 
à = lim — > log,|2z 1 
a, 82 l | 
n=1 
since a = 22 + c- Ibis not clever to do lots of 


logs, so I have reworked the formula and use 


N 
Ja ] [22 [2] 


n=1 


Then A= Inf 
N. In2 


1: you have chaos if and only if J> 1. If you are 


but you need only compare /to 
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iterating c in the center of a component, J will be 
equal to or very close to zero. In this case, evaluat- 
ing the logarithm to get A results in an error. Using 
formula [2] saves evaluating the logs in formula [1] 
and gets the same results. The factor of 2 in [2] 
tends to keep J of reasonable size during the itera- 
tion. 


How well does equation [1] or [2] work as a detec- 
tor of chaos in Mandelbrot calculations? I have 
explored this with mixed results. Along the real 
line, J/remains ‘undecided’ until n is greater than 
some value, call it & In or near a periodic compo- 
nent, I find that kis at least equal to the period. For 
example, in the period 4096 component centered 
at -1.4011551704444 Cin the Feigenbaum sequence 
approaching -1.401155189092), J— 1 remains posi- 
tive most of the time, dropping briefly below zero 
at n = 1024 and 2048. Only after n>4096 does 
J become very small. But even then it jumps above 
one a few more times, depending on how close c 
is to -1.4011551704444. I have traced the Feigen- 
baum component sequence to that with period 
262,144. Now if I were to have used the Chao- 
someter to cut off computations for chaotic points 
early, Pd have missed these small high period com- 
ponents. 


On the other hand, chaotic points on the real line 
(with A = 0.25 , for example) drive J quickly to 
very large values (J> 102% ). These can certainly 
be abandoned as chaotic saving us the trouble of 
following the iterations to some maximum iteration 
value, since the iteration will not ever diverge 

(—2 >c>-1.401155... ). For these points the test 
is useful. 


The next example I tried was c = complex points 
on the boundary of the main cardioid. These are 
either periodic (tangent to a touching bud) or cha- 
otic. 1 have posted on FRAC-L a list of chaotic c, 
some of which I tested here. The formula does not 


work as well as I would like, since lim 


N => œ 

tationally vague. Let’s start with a definitely peri- 
odic point, c = 0.25 + 0.51. After more than a 
million iterations, à approaches zero, mostly from 
the minus side, but not completely. It occasionally 
jumps positive for one or two iterations. A nearby 
chaotic point is c = 0,2850441472 + 0.46074920961. 
Again, A approaches zero, but now mostly from 


is Compu- 


the positive side, but even after more than a mil- 
lion iterations, it still drops to negative Chinting at 
periodic). 


It would appear that À for c on the edge of a com- 
ponent will approach zero in a jerky manner, and 

the test is not very useful. The test will work well 

for points on filaments, points far away from com- 
ponent edges. 


I recommend it for use in Mandelbrot calculations 
along with period detection and other tests and 
devices. Since the only filament we know well is 
the real line, the test should at least be used when 
Cis pure real. 


Implementing a 1/f Noise Algorithm 
Gabriel Landini 


Many computer simulations are based on random 
white noise generators provided in program librar- 
ies. Voss and Clarke [Vos88] demonstrated that 
there is another type of random fluctuation almost 
universally present: 1/f noise, which has a spectral 
density which is proportional to the reciprocal of 
the frequency — a random fractal noise. 


Although 1/f noise occurs commonly, there is no 
clear explanation of how it arises. Some believe 
that 1/f noise arises when processes change from 
chaotic to periodic behavior. Others attribute it to 
multiple additive or multiplicative random pro- 
cesses [Wes89]. 


Sources of 1/f fluctuation include the annual floods 
of the Nile river, sunspot activity, the fluctuations 
of electric current in thin metal sheets, cardiac elec- 
trical activity, stock exchange prices, amplitude 
fluctuations in human speech, and pitch fluctua- 
tions in music. 


In many cases computer simulations would be 
more realistic using 1/f noise rather than white 
noise. In The Science of Fractal Images, Voss gives 
a method for generating 1/f noise using the inverse 
fast Fourier transform (FFT) which requires exten- 
sive computation. 


I implemented an algorithm for 1/f noise in Quick- 


AA eee 
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BASIC which is based on an idea of Voss published 
elsewhere [Gar78] in which the values of random 
generators (dice) are added following a binary 
table of frequencies. 


1/f noise possesses a spectral density inverse to 
frequency: high frequencies, low power. This also 
means that 1/f noise possesses scaling properties: a 
1/f sequence will look statistically similar if we plot 
the same sequence skipping 1 value, 2 values, 4 
values, etc. Although this is true of white noise as 
well, it is worth trying to see what happens when a 
one-dimensional random walk or a sine wave are 
scaled. 


The program generates a sequence of length 2¢ 
depending on the number d of ‘dice’. The dice 
generate values from 1 to v, and the data have a 
Gaussian distribution of values between d and vd 


with an arithmetic mean a ue 

2 
I ran the program 1400 times with 9 dice (length = 
512) using values of v= 3, 4, 6, 8, 16, 32, and 64 
(200 times each). Each time I calculated the spec- 
tral density and its slope using the least squares 
method. The slope and correlation coefficient of 
the regression lines of a log-log plot of the spectral 
density of the 1400 runs did not show differences, 
so to get a wider range of values it is valid to 
increase v. 


Figure 1 shows the output of 1024 values of the 
algorithm analyzed with the Fast Fourier Trans- 
form. The sampling interval is one second, and 
there are eight degrees of freedom. A line with 
slope -1 — ideal 1/f noise — is superimposed for 
teterende; 


Another attractive characteristic of 1/f noise is its 
autocorrelation, which falls between uncorrelated 
white noise and the short time correlations of 
Brownian motion. 


Figure 2 shows the autocorrelation function of 1/f 
noise, which decays exponentially. 


Further experiments are possible, for example to 
create a color palette shading from white to black 
to blue to red, and then plot a sequence of pixels 
with the color value taken from the 1/f sequence. 
In this way, the autocorrelation of 1/f noise can be 
compared to that of white noise, a sine wave, or 


Brownian motion. 
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Figure 1 — Spectral Density 
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My study of clinical patterns of a buccal disease 
(periodontal disease) has some interesting theoreti- 
cal implications [Lan91a]. 


In another area, I found that 1/f data translated 
into music produces very pleasant results [Lan91b], 
as Voss discovered. 
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Quicksilver Starship 
— Rollo Silver 


Introduction 

Making a fractal movie or video of significant 
length requires a huge amount of computation. 
The successive frames of the production are usu- 
ally computed independently, although there is 
clearly a great deal of frame-to-frame redundancy. 


I present here a sketch of an algorithm — which I 
call the Quicksilver Starship algorithm, or QSA — 
which harnesses that interframe redundancy to 
provide enormous reductions in compute time. 


A Reversal 

Usually one thinks of a fractal image as a mapping 
of a region in the complex plane onto the rectan- 
gular pixel array of a monitor screen. By reversing 
this view, considering the mapping to go the other 
way— from the pixel array to the complex plane— 
we can get a significantly new insight into the pro- 
CESS: 


Concept 

Abstractly speaking, a fractal image is composed 
by computing the dwell value ACP) — a non-nega- 
tive integer — for each point Pin an mx n rectan- 
gular array of equally spaced points, and then 
coloring the corresponding pixels of an mx n 
color monitor screen according to a color mapping 
table or function T(n) which maps non-negative 
integers into colors. 
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In many cases — all those of interest to me, for 
instance — the time required to compute ACP), is 
proportional to the value ACP) itself. 


The application of OSA is not in generating a sin- 
gle fractal image, but in generating a whole series 
of closely related images displayed at high-speed 
as a series of frames to create the illusion of con- 
tinuous motion through a fractal space. 


It’s fine if we have a computer powerful enough to 
do the necessary calculations quickly enough; oth- 
erwise we have a basic problem: How can we 
make use of the information computed in one 
image of the series to reduce the complexity of 
computing the next? We need a technique which 
will bridge the gap between the brute force capac- 
ity of our computer and our requirements. 


I have gone through a series of models of decreas- 
ing complexity, arriving at last at one which is 
breathtakingly simple. All the models involve the 
idea of separating one set of points in the complex 
plane — those whose dwell is computed — from 
another set: the points corresponding to displayed 
pixels. Let’s call the former points stars and the lat- 
ter punkts. For example, the figure below shows a 
set of stars lying in one rectangular array, and a set 
of punkts (represented by disks) lying in another: 


Figure 1 — Stars and Punkts 


The “breathtakingly simple” idea is to compute the 
color value y at a punkt P using the dwell value at 
the nearest star S: y¥P) = MACS), where T and A 
are the color mapping function and dwell function 
described above. 


In order for this to work, we have to maintain a 
relationship between the two arrays. We can’t just 
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compute an array of stars (say the punkts in the 
first frame) once and for all, and then use them for 
all our punkt arrays. The problem is that there may 
not be a star sufficiently close to each punkt. 


Let's make precise what we mean by “sufficiently 

close”. Let's define the mesh tt of a punkt array to 
be the distance between adjacent punkts. We then 
require that the distance 0(P, $ between a punkt 
Pand the nearest star $ to it be no more than half 
the mesh: 


OCP, 9) <u/2 [4] 
where we define the distance not as the usual 


Euclidean distance. (X.=XJ*+ (Ve V5)? ¿ but 


as the “Hausdorttmetric”: 
OCP, 9 = max|xp— xd; | Yp- Vs) [5] 


Thus each punkt has a square neighborhood or 
box around it of “diameter” u (Figure 2, below). 


Stars whose centers are inside the box are near 
enough to the punkt; stars outside are not. 


The idea, then, is that we only need to compute a 
star dwell when there is no star sufficiently close to 
the punkt at the center of a box. To keep things 
simple, we'll keep the stars organized into a rect- 
angular array as much as possible. 


If the box around a punkt has no Star in it, we 
compute the dwell of a point within the box, turn- 
ing it into a star. We simply choose the point which 
is the center of the square. 


—— 
in ME 
SE out 
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Figure 2 — Hausdorff neighborhood of a punkt 


We're now ready for an example. 


Example: Zooming Straight In 


We start with a 7x5 array Il, of punkts, and com- 
pute the dwells for a 7x5 array 2, of stars aligned 


with II, (Figure 3). 
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Figure 3— Initial Array of Punkts and Stars 


We define H, to be the mesh i of the initial punkt 
array II). 


Let’s see what happens as we zoom in to the cen- 
tral half of the hull of this figure (Figure 4), taking 
24 frames to do so. This rate — 24 frames to halve 
the frame size — seems like a reasonable maxi- 
mum speed. 


Figure 4 — Zooming In 


Let's make the series of frame sizes a geometric 
series, so the ratio of sizes of any two successive 


frames is 74/2 = 1.0293 . Table 1 (later on) gives the 
sizes of the successive frames relative to frame 
zero, the first one. Let’s retain the stars of the origi- 
nal figure, but replace the punkts by their square 
neighborhoods (boxes) of diameter u. The initial 
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situation is shown in Frame 0. 


Frame 0— Initial Configuration 


Frame 6 shows the state of affairs after six zoom-in 
steps. The punkt array has shrunk to 0.8409 of its 
original size: the grey circles cover stars that have 
just passed outside the frame, and the circled stars 
have just been created at the centers of empty 
boxes. 


Frame 6b shows the frame after the escaped stars 
and the circles around the new stars have been 
removed. 


Frame 8 shows the state of affairs after zoom step 
8. The punkt array has shrunk to 0.7937 of its orig- 
inal size. The grey circles cover stars which have 
moved outside all the boxes. In addition the 16 cir- 
cled stars have been created to fill empty boxes. 
Frame 8b shows frame 8 as it actually appears. 


Frame 8 — After eight steps, with escaped stars 
and new stars 


Frame 8b — stars and boxes 


Frame 6 — After Six Steps, showing escaped stars 
and new stars 


A na A 
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Frame 12 — After 12 steps, with escaped stars and 
new stars 


Frames 12 and 12b show the state of affairs after 
zoom step 12. The punkt array has shrunk to 
0.7071 of its original size. Ten stars have moved 
out of the frame, and four stars have been added at 
the centers of empty squares. 


Frame 12b — stars and boxes 


Frames 14 and 14b show the state of affairs after 
zoom step 14. The box array has shrunk to 0.6674 
of its original size. Six stars have been added at the 
centers of empty boxes. 


Frame 14 with new stars 
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Frame 14b — Stars and boxes 


The stars which appear to be on the boundary 
between two boxes, or in the common vertex of 
four boxes, are actually not: the arrows indicate 
which box such a star is really in. 


Instead of continuing to show the details of stars 
lost outside the shrinking box array and new stars 
created to occupy newly empty boxes for the suc- 
cessive frames, we summarize the situation in the 
following table, and then show the final frame: 
#24. The table has entries for only those frames 
which have significant activity, i.e. which have 
either lost stars or have gotten new stars created in 
them. 


Table 1 — Frames with significant activity 


Scale Lost New 
Factor stars stars 


0 1.000 0 35 
0.8409 10 10 
8 0.7937 10 16 
12 0.7071 10 4 
14 0.6674 0 6 
15 0.6484 0 6 
16 0.6300 10 8 
17 0.6120 10 8 
18 0.5946 4 0 
20 0.5612 6 4 | 
21 0.5453 4 0 | 
23 0.5147 0 4 


0.5000 
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Frame 24 — lost stars and new stars 


Let’s sum up. If we did not take advantage of the 
frame-to-frame redundancy of the process, we 
would have computed 35 stars for each of the 24 
frames of the zoom shot after the first: a total of 
840 stars. This requires a computer having star- 
compute power of 35 stars/frame. 


Applying OSA, we need only compute 74 stars for 
the 24 frames, requiring star-compute power of 
3.08 star/frame: an elevenfold reduction. 


Practical Situations 

In practice we shouldn't take the speedup factor of 
eleven too seriously: it will be affected by the size 
of the punkt array, the zoom rate, and panning. 


Pd expect the speedup factor to increase with 
increasing array size: In proportion to the linear 
dimension of the array, I would suggest, since for 
an mx n array, the number of stars is proportional 
to mn, while the number of stars lost per frame 
should be proportional to the number of boundary 
pixels: 2(m+n-—A4) . 


Pd expect the speedup factor to decrease with 
increasing zoom rate. If we define the “half-life” 7 
to be the number of frames computed between 
one frame and the following frame half its size 

(1 24 for the-examole-eiven), J dvexpect the 
speedup factor to be proportional to T. 


Pd expect the speedup factor to decrease with 
increasing pan rate. If we define the pan-time to be 
the number of frames computed during the time it 
takes for a star to move from one edge of the 
punkt-rectangle to the opposite edge, Pd expect 
the speedup factor to be proportional to the pan- 
time. 


For example, if we produce a fractal zoom movie 
which runs at 24 frames/second, and we have a 
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rectangle of 640x480 pixels, and we zoom in 
toward the frame center at a rate which halves the 
frame-size every two seconds, I would expect a 
speedup factor of about 2,000! 


If such staggering speedup factors prove in fact to 
be the case, we will find ourselves at the dawn of a 
new era, in which we can produce not only pre- 
computed, canned displays on a personal com- 
puter, but also undertake real-time journeys 
through fractal space, controlled by a human oper- 
ator using a joystick and throttle! 


Implementation 


If you are interested in implementing QSA, you 
should contact me if you need help with details of 
the algorithm or data structures to support it. In 
any case, we should discuss the conditions of your 
implementation (e.g. freeware, shareware, com- 
mercial application), and arrange for permission 
for you to use my algorithm. 


Complex Inverse Trig Functions 
Kerry Mitchell 


On Compuserve, a question came up about inverse 
trig functions of complex variables. These formu- 
las, from Complex Variables and Applications by 
R.V. Churchill and J.W. Brown (McGraw-Hill, 
1984), may be useful to die-hard fractal freaks such 
as myself. 


asinz = —1log (t{z+ y1-— 27) 
acosz = —1log (z+i1— 22) 


atan Z 


II 
NI s. 
O 

QQ 


asinh z = log (z+A/2?+1) 
acoshz = log (z+ J22- 1) 


i l+z 
2:98 


ll Z 


atanh z 
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RS adds: 


For purposes of calculation, these formulas must 
be reduced to functions of a real variable and sim- 
ple complex arithmetic, so we need the following 
formulas as well: 


Jz = let cose tale . 
2 2 

logz = logr+1(0+279) , 
ME gag ley de Lo 3e3 


Where 7S Wee. Dos atan, and of course 
ZA AA 


Note that the square root and log functions are 


multi-valued: Jz has two values, while logz has 
an infinity of values, one for each integer value of 
n. Which of the multi-values to use when calculat- 
ing is a deep question, beyond the scope of this 
note. For some insight into this conundrum see the 
article Bifurcations, in Amygdala #3. 

For those readers who may have some difficulty it 


seeing how to turn the math into code, here’s how 
the first function might be implemented: 


asin z = —tlog (z+ yl — 2?) 


logz = logr+1(0+2T MH) 


Jz = +r sa omen 
| 2 2 


Here's code to compute u and v given x and y, 
where “+ iv = asin (x+ iy) -TIl assume nto be 
externally defined for use in the log function. pt is 
te constant i=" 3. ss 

xl = 1- x42 + y*2 

VL = =Z*X%*Y 

GS Sart sarta a yL) 

E = atan (is xy 72 

x2 = sy + q*cos(t) 


Vo. = xe Ge sac) 
atan ly2/x2) + 2*pisn 
=109 (202 + Y2N 2) L2 


¡e 
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Book and CD Reviews 


“Fractal Fantasies” by Lee H. Skinner 
(A CD ROM) 


A review by John Jones 


This is not the first time I’ve seen fractal images, 
but it is the first time I’ve used a CD-ROM. The first 
thing I noticed was that a CD-ROM holds a lot of 
data — there are over 2,000 images on the disk. 
And each image actually appears three times — in 
1024-by-768 resolution, 640-by-480 resolution, and 
16-to-a-screen thumbnail resolution. I viewed the 
ROM on a PC with a VGA video card, so I looked 
at the 640-by-480 images, then transferred the 
high-resolution versions of a few of my favorites to 
a SparcStation for detailed examination. 


But who would want 2,000+ fractal images on CD- 
ROM? Only a hard-core fractal enthusiast. And such 
enthusiasts often prefer to create their own images. 


Fortunately, the CD also includes source code for 
FractInt. So it doesn’t just offer a passive viewing 
experience; the images provide a guide to the 
enormous range of possibilities in FractInt itself, 
and may inspire the viewer to explore indepen- 
dently. 


The images don’t provide a map to the worlds 
reachable via FractInt; rather, they are snapshots 
from one person’s extensive wanderings through 
those worlds. And on the whole they are snap- 
shots, not Ansel Adams masterpieces. There are 
masterpieces on the disk: I have a list of forty 
images that are as impressive as anything I’ve ever 
seen. EPSCR/1ASPENEP.GIF, for example, is a beau- 
tifully composed picture, which reminds me of 
Escher's Three Worlds, RCL10/1THEBLOB has a mar- 
velous color scheme, showing fluid orange 
pseudopods reaching up under a canopy of fir 
trees. And there are several hundred others which 
are sufficiently interesting to make me want to 
explore the regions that they came from. But per- 
haps 80% of the images left me cold. This is partly 
a matter of familiarity — there are whole series of 
images that have the same texture as the Mandel- 
brot set, and I found these far less interesting than 
textures new to me, such as the Kandinskian com- 
positions in the E1E2 subdirectory — and partly a 
matter of the care that went into crafting each 
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image — some of the pictures don't seem to have 
a focus of interest, while others have strangely 
clashing colors. ZEPE-F/1ZEXPF90.GIF, for example, 
is pleasing but poorly framed, while PARABLA.GIF is 
marred by the contrast between juxtaposed areas 
of blue and orange. 


This does not detract from the usefulness of the 
CD; sometimes a snapshot is enough. And other 
viewers may find gems among the 80% I’d dismiss. 


My most serious complaint about the CD is the dif- 
ficulty of finding where a particular image came 
from. Sometimes this is not too difficult — for 
example, the images in directory CHCO are 
described as of type FZPPCHCO, and searching the 
FRACTINT.FRM file of the PROGRAMS subdirectory, I 
was able to find the corresponding formula. But 
the images for which I most wanted the source, the 
beautiful and mysterious images of MBCOMP, had 
only the tantalizing description ‘Mandelbrot Com- 
parison Images’, and I still have no idea how to re- 
create them. Maybe this information is stored in 
FractInt. But many potential purchasers of the CD 
may have a computer, like my Sparc, on which 
FractInt won't compile. So if the slide-show is to 
stand alone, it would be nice to have this informa- 
tion ready to hand. A hypertext interface might be 
the best way of doing this. 

A second flaw in the CD is the fact that several 
subdirectories are missing from the 640 image 
series. In particular, HALLEY, LAMB, LEES3D and LPC 
are featured at 1024 and in thumbnail resolution, 
but not at 640. This is a pity, because the 3-D pic- 
WES are Excelent: 

My final peeve is that ‘ancestor’ is misspelled as 
‘ancester’ throughout. 

In summary, I would say that this CD is certainly 
worth its price. Viewing it has suggested many ave- 
nues for further exploration. 

— John Jones 


$35.55 + $5 shipping, from Amyedala. 


“Fractals and Chaos in Geology and 
Geophysics” 
A review by lan Lumb 


Fractals and Chaos in Geology and Geophysics, 
Donald L. Turcotte, Cambridge University Press, 
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1992 pp x+221. ISBN 0-521-41270-6; 
QE33.2.M3T87. Price: $54.95 (hard cover). 


This book is intended to introduce the reader to 
the concepts of fractals and chaos in geological 
and geophysical contexts. Designed to serve the 
senior-undergraduate level, this text is supple- 
mented by an up-to-date reference list plus glossa- 
ries of terms, units and symbols; 15 of the 17 
chapters include instructive problems, many of 
which have solutions listed at the end of the text. 
Although the text is ideally suited for those well 
versed in geology and geophysics, it is sufficiently 
self-contained to be accessible to those from other 
disciplines. 


After a brief introductory chapter, the next seven 
chapters focus on fractals. In the second chapter, 
the concept of a fractional dimension is clearly 
introduced in one, two and three dimensions. The 
author moves swiftly from a consideration of the 
fractal dimension of ‘classical’ fractals, to ‘real- 
world’ coastline and topographic data. Although 
the use of power-law Statistics in geology and geo- 
physics predates the fractal concept (B. B. Mandel- 
brot, 1967), it is instructive to explore a variety of 
size-frequency distributions using a fragmentation 
framework provided by fractals. Since fragmenta- 
tion involves the highly nonlinear process of frac- 
ture propagation, and interactions between these 
fractures of varying scales, a fractal-based descrip- 
tion seems to be called for; the models of fragmen- 
tation developed in Chapter 3 allow for a 
discussion of a variety of processes which include 
explosions/impacts, weathering, and porosity. 
Starting with the classic Gutenberg-Richter law, 
fragmentation models are used to illustrate the 
scale invariance of frequency-size distributions of 
earthquakes in Chapter 4; the same scale indepen- 
dence of volcanic eruptions closes out this chapter. 
In the next chapter, the relationship between mean 
grade and cumulative ore tonnage is illustrated for 
the production of several economically interesting 
ores, as well as for oil field production. When frac- 
tal distributions are related to probability, one can 
study the clustering of events in time or space; this 
is the subject of Chapter 6, where a brief applica- 
tion to seismicity data is considered. 


Up to this point, the author indicates that only self- 
similar fractals have been considered; i.e. those 
that are statistically similar at different scales, save 
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for a multiplicative factor. On the other hand, self- 
affine fractals involve scaling factors raised to a 
power, and are usually treated with spectral meth- 
ods. Thus in Chapter 7, the author presents a com- 
prehensive review of fractal-dimension estimation 
via spectral methods in one and two dimensions; 
real-word examples include seismic reflection coef- 
ficient, topography, and geoid data. Four color 
plates close off this chapter with illustrations of 
synthetic topography, plus fractal and roughness 
parameter characterization of digitized topography. 
The scale invariance of erosive processes allows 
for landscape evolution modeling, via spectral 
methods, in Chapter 8. 


In Chapter 9, an introduction to dynamical systems 
is provided through an analysis of the van der Pol 
equation; a solid review of the ‘classic’ logistic map 
follows in Chapter 10. The author illustrates how 
the logistic map iteratively evolves as based on dif- 
ferent choices of a constant parameter and an ini- 
tial value; in so doing, fixed points, limit cycles, 
chaotic behavior, and bifurcations are all well illus- 
trated. Quantification of the degree of chaos, as 
provided by the Lyapounov exponent, is also 
explained and illustrated in this same context. 


The production of an earthquake involves the for- 
mation of a pressure bond between adjacent 
crustal blocks, which allows for a build up of elas- 
tic strain potential energy. Once the potential 
energy exceeds the strength of the pressure bond, 
brittle fracture occurs to release the strain, and the 
potential energy is converted into kinetic energy. 
In Chapter 11, a slider-block model is used to 
describe the kinematics of this earthquake process. 
This model consists of a block of mass M, which is 
driven by a constant-velocity driver at speed V 
until the force, kx, on a spring exceeds the static 
friction force. An equation of motion, which 
applies once sliding has begun, is used to describe 
the time evolution of the system. Single-block solu- 
tions provide the foundation for the more ‘Earth- 
like’ wwo-block analyses that are also discussed in 
the text. It is in these latter analyses that the signs 
of chaos can be well illustrated, and promising 
comparisons with real earthquake data can be 
made. 


In his classic article (Deterministic Nonperiodic 
Flow, J. Atmos Sci., 20, 130-141, 1963) E. N. Lorenz 
essentially ‘discovered’ what we now refer to as 
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chaos, i.e. an exponential sensitivity to initial con- 
ditions. Lorenz was interested in the study of ther- 
mal convection in a fluid layer heated from below. 
In Chapter 12 a thorough review of Lorenz’s prob- 
lem, which involves a set of three coupled total 
differential equations, is provided. The author pro- 
vides the details for the derivation of the Lorenz 
system, the calculation of the system's evolution, 
and projections of the resulting Lorenz attractor. In 
the following chapter, the author uses the formal- 
ism developed for the Lorenz problem to investi- 
gate the nature of thermal convection in the Earth's 
mantle. Using a double Fourier series expansion 
and scalings appropriate to mantle convection, the 
author illustrates the chaotic nature of this hydro- 
dynamic instability. In Chapter 14, a numerical 
solution of the Rikitake two-disk dynamo equa- 
tions is shown to have some properties in common 
with the solution of the Lorenz equations, and pos- 
sess some similarities to actual sequences of mag- 
netic field reversals recorded over geologic time 
scales. 


The last two chapters, which precede the very brief 
concluding chapter, are concerned with large inter- 
active systems. In Chapter 15, the renormalization 
group method is applied to fluid flow in a porous 
medium, fragmentation, and a fractal tree model 
for rupture of a fault. Self-organized criticality is 
used to re-investigate the earthquake frequency- 
size distribution problem, and seems to show some 
predictive promise. 


This is a purposefully, concisely, and well-written 
book on an important formalism for geology and 
geophysics. It is fortunate that someone with the 
author's experience and stature has taken the time 
to put together this impressive synthesis. Despite 
its relatively small length, this book spans many 
interesting topics and is completely self-contained. 
It would thus be highly suitable as a senior under- 
graduate or graduate-level text. Although the over- 
all organization of the text is satisfactory, chapter 
subheadings would be a welcome addition. 


Those interested in this subject might also want to 
examine the monograph Nonlinear Dynamics and 
Predictability of Geophysical Phenomena, Geo- 
physical Monograph 83, American Geophysical 
Union, 1994, which is co-edited by this same 
author. 
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lan Lumb, Research Scientist 

KelRescarch Corporation 

Consultants in Atmospheric Science & Technology 
850-A Alness Street, Suite 9 

Downsview, Ontario M3] 2H5 CANADA 

Voice: 416-736-0521; FAX: 416-661-7171 

Internet: <kel@nexus.yorku.ca> 


“Chaos In Wonderland”, A Book By 
Clifford Pickover 


From the publisher’s release: 


“Alice in Wonderland meets 2001: A Space Odys- 
sey! 


“Stir together a mixture of fractals, chaos, computer 
eraphics, and science fiction. What do you get? 
You get a dazzling introduction to chaos science 
by Clifford Pickover, IBM's indefatigable computer 
scientist. Dr. Pickover alleges that his gorgeous 
swirling art, generated by the strange attractors of 
two simple formulas are the dreams of limbless, 
brainy creatures living below the icy surface of 
Ganymede, Jupiters largest moon. 


“Things flow about here! exclaimed Alice, in the 
queer little shop behind the mirror; but never in 
Lewis Carroll’s wildest imagination could she have 
dreamed of the flowing patters in Pickover’s Won- 
derland.” 

—Martin Gardner, Scientific American 


“A fascinating project, a worthy addition to Flat- 
land and The Planiverse? 


Arthur C. Clarke 


“Pickover weaves a fascinating and entertaining 
tale of fact and fiction in Chaos in Wonderland. 
While learning of the intricacies of the world of the 
Latocoecarfians, we are skilfully introduced to the 
mathematics of chaos. His mixture of fantasy and 
mathematics creates an ceriness of reality that is 
both informative and thought provoking.” 
—Theoni Pappas, author of The Joy of Mathematics 


“A leading expert in computer visualization, Pick- 
over now enthralls us with his art, mathematical 
games, and science fiction. This latter-day Lewis 
Carroll introduces us to alien creatures with com- 
puter brains and mathematical souls. Their social 
status is based on the beauty of geometrical pat- 
terns communicated to one another with infrared 


beams. You will be delighted at the mathematical 
wonderland Pickover provides using media more 
suited to us humans. 


—Prof. J.C. Sprott, author of Strange Attractors 


“Chaos Under Control” 
—A.G. Davis Philip 


I have before me a recently published book, 
“Chaos Under Control, The Art and Science of 
Complexity” by David Peak and Michael Frame of 
Union College. A short summary of the contents 
can be given by listing the subjects of the ten chap- 
LETS: 

1. Introduction — why chaos and fractals are inter- 
esting 

2. Iterated Function Systems 

3. Fractal Dimensions 

4, Chaotic Dynamics and Graphical Iteration 

5. Order in Chaotic Dynamics —Bifurcation Dia- 
grams 

6. Graphical Test of Chaos 

7. Mandelbrot set and Julia sets 

8. Fractal Basin Boundaries 

9. Cellular Automata, Neural Nets, Sandpiles and 
Artificial Life 

10. Fractals and Chaos in the Humanities and 
Social Sciences. 


There are 16 color plates in the center of the 
book., which is now being used in courses at 
Union College, Yale University and the University 
of Richmond. In a foreword by Mandelbrot, he 
writes, “The observation is that fractals — together 
with chaos, easy graphics and the computer — 
enchant many young people, which — in turn — 
makes them excited about learning mathematics 
and physics.... In this book David Peak and 
Michael Frame have taken a step in realizing this 
hope.” 


“Fractal Cities” 
—Michael Batty and Paul Longley 


From the publishers description: 


Fractal Cities is a study of the development and use 
of fractal geometry for understanding and planning 
the physical form of cities, showing how this 
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geometry enables cities to be simulated through 
computer graphics. It shows how cities evolve and 
grow in ways that at first sight appear irregular, but 
which, when understood in terms of fractals, illus- 
trate an underlying order that reveals their com- 
plexity and diversity. 


The book contains sixteen pages of computer 
graphics and explanations of how to construct 
them, as well as new insights into the complexity 
of social systems. The authors provide an intelligi- 
ble and gentle introduction to fractal geometry as 
well as an exciting visual understanding of the 
form of cities, thus providing one of the best intro- 
ductions to fractal geometry available for non- 
mathematicians and social scientists. 


Fractal Cities can be used as a text for courses on 
geographic information systems, urban geogra- 
phy, regional science and fractal geometry. Plan- 
ners and architects will also find that there are 
many aspects of fractal geometry in this book rele- 
vant to their own interests. Furthermore, those 
involved in fractals and chaos, computer graphics 
and systems theory will find important methods 
and examples which are germane to their work. 


Michael Batty is Director of the National Center for 
Geographic Information and Analysis in the State 
University of New York at Buffalo, and Paul Long- 
ley is Reader in Geography at the University of 
Bristol. 


Book Reviews Wanted 

If you run across a book on the subject of fractals, 
etc. that seems particularly illuminating, and has 
not yet been reviewed in Amygdala, think about 
writing a review of the book and sending it in for 
possible publication in Amygdala. 
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Quaternion Square Root 
Rollo Silver 


This note on quaternion arithmetic arose in 
response to a question posed by Bill Bell on the 
Internet FRAC-L forum. He wrote, “I'm trying to 
write a program to display 3-D Julia sets using 
quaternions. My problem is, the math still baffles 
me a little bit, in particular trying to figure out how 
to calculate the square root of a quaternion. I've 
looked up most of the magazine articles treating 
this subject, but their explanations are still not pen- 
etrating my thick skull. I keep winding up trying to 
find the arctangent of a complex number, which 
I'm not sure is possible.” 


I worked out a simple algorithm for the square 
root of a quaternion, not involving anything as 
complicated as the arctangent of a complex num- 
ber. Gerald Edgar pointed out errors in my first 
two attempts to get the algorithm right— Now I 
think it’s right! 


So here we go, computing the quaternion r which 


is the square root r = Ja of a quaternion 4. 


Let 

q = Qt 14, + q, + RQ; [6] 
and let 

r= +17, +31, + kr, [7] 
bea TOOL OF 

ye => [8] 


In summary we have: 

(1) If q = 0, then [8] has precisely one root: 

em ie 

(2) If q is negative real, then [8] has an infinity of 


roots [7] with n= O and 77+ 7+ 7 = =o 


(3) Otherwise [8] has precisely two roots, 
E sta kq; 


2p 2p 2p 


dot (G+ G+ G+ B 
(r. 


Now for the gory details. 


) -where 
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From [7] we have 
P= (+1, +]1,+ kr) (+11, +Jmt+kn) , 
which expands and simplifies to 


A E Be ye 
y = re — v2 — 18 — 2 + DIN, 1, + 2/1,1, + 2 kr,7, [9] 


Substituting [9] and [6] in [8], we get 
R—R-%-R+21T4 +21, + 2RN; = 
4 + iq, + jq, + RG, . Equating the four quaternion 
components separately, we get: 
i—-A-R-%=4 (101 
and 
27% = Gs 2%% = I, 277% = q [11] 


Our task is to find all real iTo %, % which satisfy 
[10] and [11]. 


Case 1: q is real. 
That is, q = qa. We have two subcases, % = 0 


and 7 #0. 


Case 1.1: q is real (4 = q,) and 7, = 0. 
From [10], 77+ 722+ 1% = -q,andso qa £0, and 
any 1, %, Y, Satisfying 77+ 13+ 1% = =q define a 
root r = ir, +Jm+ kr, of [8]. For example, if 

q = —1, the following are among the infinity of 
+ ¡+ k 


— A 


roots of [8]: 47, R, 


Case 1.2: q is real (4 = q, ) and r7,+0. 
Since % #0, while each of q,, q», q3 is 0 (since q 
is real), [11] implies that r, = n = % = 0, in 
which case [8] reduces to 7 = q,. We must have 
qd) 20, and [8] has either one or two roots, 
depending on whether q, = 0 or q,>0.If q = 0 
then r= 0. 1f g>O then r= +q. 


Case 2: q is not (pure) real; i.e. not all 
of g,, q, gq, are zero. 


In this case % #0, lest by [11], q, = 0, q, = 0, 


q = 0. We can then divide each of the equations 


[11] by r, and get 


_ A _ 4) _ Ww 
eal ee ee [12] 


Substituting [12] in [10]: 


n-(#)-(#)-(#) - a 
AD) 2H) eG 
4% -4q01- (gq? + G+ q3) = 0, whose roots are 


Qt J+ G+ Gt a 


— 
nonzero, 7>0. But g< ql + q? tatg? so we 


R = . Since 7 is real and 


must take the sign of the square root to be posi- 
tive: 


Pe do+ NG + ++ 113] 
A ae cali 


We have p and —p as solutions of [13], where 


Got AGB + a? + G+ GF 


Applying [7] and [12], in this case [8] has precisely 
two roots: 


id Ib, = 


2p * 2p * ap [15] 


r= a(p+ 
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ham University on several applications of fractal 
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Among his other interests are math games, chess, 
Go, and hiking. 


lan Lumb’s research efforts are geared towards the 
completion of his Ph.D, degree in Physics:and 
Astronomy at York University (North York, 
Ontario, Canada), and towards several contracts at 
the KelResearch Corporation (Downsview, 
Ontario, Canada). It was at KelResearch that he 
first had the opportunity to delve into the depths 
chaos and fractals, and continues to study image 
processing and analysis via fractal-based methods. 
lan also teaches at York University and is an avid 
Internetworker. He is married with three daugh- 
ters. 


Each One Bring One 

Subscribers, I need your help! The more Amygdala 
subscribers we have, the better the newsletter will 
be! 


I have been working to build circulation by post- 
ing notices on Internet forums and newsgroups, 
and on Compuserve forums, but I need your help! 
Please reach out among your friends, acquaintan- 
ces, coworkers, etc. to sign up new subscribers for 
Amygdala. 


The subscription rates for ten issues are $28.20 
(US), $32.00 (Canada & Mexico), $39.00 (Other). 
First Class or Air Mail postage included. 


We also need material for future issues! My writing 
takes the lion’s share of this issue, and I think that 
we need more diversity. If you can write, and have 
material that you think would interest Amygdala 
readers, please submit it for publication. See the 
Submitting Articles guidelines on this page. 


Amygdala is published by Rollo Silver every 
eight to ten weeks. 

Address articles, letters, comments, subscription 
orders, to: Amygdala; Box 219; San Cristobal, 
NM 87504; USA. 

Email: rsilverObeta.lanl.gov 

Phone/VAX: 505/586-0197 

Compuserve: 71174,1453 


This issue was produced using FrameMaker on a 
Macintosh Mex computer. 


DA Amygdala #38 


Submitting Articles 


Here are some guidelines for submitting articles for 
publication in Amygdala. 


1. Type your article as you would like to see it 
appear in the newsletter. Please do not send hand- 
written drafts! 


Send it (my order of preference): (A) On 3.5” dis- 
kette in Macintosh format for FrameMaker, 
MacWrite, WriteNow, or MS Word; or as a text file. 
Please also enclose paper copy so I can see your 
intent. (B) On 3.5” diskette in IBM format, text file. 
I have no way to deal with 5.25” diskettes. (C) 
Paper copy. 


2. Illustrations should be either greyscale (suitable 
for halftoning) or black/white; not color! CA) Nor- 
mally, illustrations will be printed in full column 
width, so you should make them 3.25” wide, if 
possible — provided that they’re 300 dpi resolu- 
tion. If they’re grainier, make them larger if possi- 
ble, so that they’ll look good when reduced to 300 
dpi. (B) Make sure that you clearly indicate which 
illustrations go where in the text! (C) All in all, it’s 
better not to have captions welded into your pic- 
tures. Let me put them in ad lib. (D) I can handle 
illustrations on diskette in MacPaint, MacDraw, 
CricketDraw, Photoshop, or Adobe Illustrator for- 
mats. 


3. Please send along a short biographical note, 
which I will try to publish in the same issue as 
your article. 


4. Please include your telephone number, in case 
I have to reach you in a hurry with questions. 


Circulation/Renewal 


As of October 29, 1994, Amygdala has 438 sub- 
scribers. 


The address used for mailing this issue to you has 
a small notation at the lower right indicating when 
your subscription expires. 


A33 or A34 indicates that this is your last issue. 
A35 orA36 indicates that your subscription expires 
with the next issue. In either case, please use the 
form on the other side of the address page to 
renew. 
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